ABSTRACT
In 2019, novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) began infecting humans, resulting in the COVID-19 pandemic. While the push for development of vaccines has yielded some positive results, the emergence of additional variants has led to concerns surrounding sustained vaccine effectiveness as the variants become the dominant strains. This work was undertaken to develop peptide-based antivirals capable of targeting both the wildtype (WT) heptad repeat 1 (HR1) domain of SARS-CoV-2 and the new HR1 variants which have developed. In silico protein mutagenesis, structural characterization, and protein–protein molecular docking were utilized to determine molecular interactions which facilitated binding of peptide-based antivirals targeting the HR1 domains. Molecular dynamics simulations were utilized to predict the final binding affinities of the top five peptide inhibitors designed. This work demonstrated the importance of hydrophobic interactions in the hydrophobic gorge and in the rim of the HR1 domain. Additionally, the placement of charged residues was shown to be essential in maximizing electrostatic interactions. The top five designed peptide inhibitors were all demonstrated to maintain good binding affinity to the WT and the variant HR1 SARS-CoV-2 domains. Therefore, the peptide inhibitors designed in this work could serve as potent antivirals which are effective in targeting both the original SARS-CoV-2 and the HR1 variants that have developed.
ABSTRACT
Coronaviruses are viral infections that have a significant ability to impact human health. Coronaviruses have produced two pandemics and one epidemic in the last two decades. The current pandemic has created a worldwide catastrophe threatening the lives of over 15 million as of July 2020. Current research efforts have been focused on producing a vaccine or repurposing current drug compounds to develop a therapeutic. There is, however, a need to study the active site preferences of relevant targets, such as the SARS-CoV-2 main protease (SARS-CoV-2 Mpro), to determine ways to optimize these drug compounds. The ensemble docking and characterization work described in this article demonstrates the multifaceted features of the SARS-CoV-2 Mpro active site, molecular guidelines to improving binding affinity, and ultimately the optimization of drug candidates. A total of 220 compounds were docked into both the 5R7Z and 6LU7 SARS-CoV-2 Mpro crystal structures. Several key preferences for strong binding to the four subsites (S1, S1', S2, and S4) were identified, such as accessing hydrogen binding hotspots, hydrophobic patches, and utilization of primarily aliphatic instead of aromatic substituents. After optimization efforts using the design guidelines developed from the molecular docking studies, the average docking score of the parent compounds was improved by 6.59 -log10(Kd) in binding affinity which represents an increase of greater than six orders of magnitude. Using the optimization guidelines, the SARS-CoV-2 Mpro inhibitor cinanserin was optimized resulting in an increase in binding affinity of 4.59 -log10(Kd) and increased protease inhibitor bioactivity. The results of molecular dynamic (MD) simulation of cinanserin-optimized compounds CM02, CM06, and CM07 revealed that CM02 and CM06 fit well into the active site of SARS-CoV-2 Mpro [Protein Data Bank (PDB) accession number 6LU7] and formed strong and stable interactions with the key residues, Ser-144, His-163, and Glu-166. The enhanced binding affinity produced demonstrates the utility of the design guidelines described. The work described herein will assist scientists in developing potent COVID-19 antivirals.